tutorial surface emg detection in space and time best practices

Abstract This tutorial is aimed to non-engineers using, or planning to use, surface electromyography (sEMG) as an assessment tool in the prevention, monitoring and rehabilitation fields. Its first purpose is to address the issues related to the origin and nature of the signal and to its detection (electrode size, distance, location) by one-dimensional (bipolar and linear arrays) and two-dimensional (grids) electrode systems while avoiding advanced mathematical, physical or physiological issues. Its second purpose is to outline best practices and provide general guidelines for proper signal detection. Issues related to the electrode-skin interface, signal conditioning and interpretation will be discussed in subsequent tutorials

Specificity of surface EMG recordings for gastrocnemius during upright standing.

The relatively large pick-up volume of surface electrodes has for long motivated the concern that muscles other than that of interest may contribute to surface electromyograms (EMGs). Recent findings suggest however the pick-up volume of surface electrodes may be smaller than previously appreciated, possibly leading to the detection of surface EMGs insensitive to muscle activity. Here we combined surface and intramuscular recordings to investigate how comparably action potentials from gastrocnemius and soleus are represented in surface EMGs detected with different inter-electrode distances. We computed the firing instants of motor units identified from intramuscular EMGs detected from gastrocnemius and soleus while five participants stood upright. We used these instants to trigger and average surface EMGs detected from multiple skin regions along gastrocnemius. Results from 66 motor units (whereof 31 from gastrocnemius) revealed the surface-recorded amplitude of soleus action potentials was 6% of that of gastrocnemius and did not decrease for inter-electrode distances smaller than 4 cm. Gastrocnemius action potentials were more likely detected for greater inter-electrode distances and their amplitude increased steeply up to 5 cm inter-electrode distance. These results suggest that reducing inter-electrode distance excessively may result in the detection of surface EMGs insensitive to gastrocnemius activity without substantial attenuation of soleus crosstalk

Recruitment of motor units in the medial gastrocnemius muscle during human quiet standing: is recruitment intermittent? What triggers recruitment?

The recruitment and the rate of discharge of motor units are determinants of muscle force. Within a motoneuron pool, recruitment and rate coding of individual motor units might be controlled independently, depending on the circumstances. In this study, we tested whether, during human quiet standing, the force of the medial gastrocnemius (MG) muscle is predominantly controlled by recruitment or rate coding. If MG control during standing was mainly due to recruitment, then we further asked what the trigger mechanism is. Is it determined internally, or is it related to body kinematics? While seven healthy subjects stood quietly, intramuscular electromyograms were recorded from the MG muscle with three pairs of wire electrodes. The number of active motor units and their mean discharge rate were compared for different sway velocities and positions. Motor unit discharges occurred more frequently when the body swayed faster and forward (Pearson R = 0.63; P < 0.0001). This higher likelihood of observing motor unit potentials was explained chiefly by the recruitment of additional units. During forward body shifts, the median number of units detected increased from 3 to 11 (P < 0.0001), whereas the discharge rate changed from 8 ± 1.1 (mean ± SD) to 10 ± 0.9 pulses/s (P = 0.001). Strikingly, motor units did not discharge continuously throughout standing. They were recruited within individual, forward sways and intermittently, with a modal rate of two recruitments per second. This modal rate is consistent with previous circumstantial evidence relating the control of standing to an intrinsic, higher level planning proces

Online tracking of the phase difference between neural drives to antagonist muscle pairs in essential tremor patients

Transcutaneous electrical stimulation has been applied in tremor suppression applications. Out-of-phase stimulation strategies applied above or below motor threshold result in a significant attenuation of pathological tremor. For stimulation to be properly timed, the varying phase relationship between agonist-antagonist muscle activity during tremor needs to be accurately estimated in real-time. Here we propose an online tremor phase and frequency tracking technique for the customized control of electrical stimulation, based on a phase-locked loop (PLL) system applied to the estimated neural drive to muscles. Surface electromyography signals were recorded from the wrist extensor and flexor muscle groups of 13 essential tremor patients during postural tremor. The EMG signals were pre-processed and decomposed online and offline via the convolution kernel compensation algorithm to discriminate motor unit spike trains. The summation of motor unit spike trains detected for each muscle was bandpass filtered between 3 to 10 Hz to isolate the tremor related components of the neural drive to muscles. The estimated tremorogenic neural drive was used as input to a PLL that tracked the phase differences between the two muscle groups. The online estimated phase difference was compared with the phase calculated offline using a Hilbert Transform as a ground truth. The results showed a rate of agreement of 0.88 ± 0.22 between offline and online EMG decomposition. The PLL tracked the phase difference of tremor signals in real-time with an average correlation of 0.86 ± 0.16 with the ground truth (average error of 6.40° ± 3.49°). Finally, the online decomposition and phase estimation components were integrated with an electrical stimulator and applied in closed-loop on one patient, to representatively demonstrate the working principle of the full tremor suppression system. The results of this study support the feasibility of real-time estimation of the phase of tremorogenic neural drive to muscles, providing a methodology for future tremor-suppression neuroprostheses

Biomechanical analysis of body movements of myoelectric prosthesis users during standardized clinical tests

Objective: The objective clinical evaluation of user's capabilities to handle their prosthesis is done using various tests which primarily focus on the task completion speed and do not explicitly account for the potential presence of compensatory motions. Given that the excessive body compensation is a common indicator of inadequate prosthesis control, tests which include subjective observations on the quality of performed motions have been introduced. However, these metrics are then influenced by the examiner's opinions, skills, and training making them harder to standardize across patient pools and compare across different prosthetic technologies. Here we aim to objectively quantify the severity of body compensations present in myoelectric prosthetic hand users and evaluate the extent to which traditional objective clinical scores are still able to capture them. Methods: We have instructed 9 below-elbow prosthesis users and 9 able-bodied participants to complete three established objective clinical tests: Box-and-Blocks-Test, Clothespin-Relocation-Test, and Southampton-Hand-Assessment-Procedure. During all tests, upper-body kinematics has been recorded. Results: While the analysis showed that there are some correlations between the achieved clinical scores and the individual body segment travel distances and average speeds, there were only weak correlations between the clinical scores and the observed ranges of motion. At the same time, the compensations were observed in all prosthesis users and, for the most part, they were substantial across the tests. Conclusion: The sole reliance on the currently available objective clinical assessment methods seems inadequate as the compensatory movements are prominent in prosthesis users and yet not sufficiently accounted for

Intramuscular EMG-driven musculoskeletal modelling: towards implanted muscle interfacing in spinal cord injury patients

OBJECTIVE: Surface EMG-driven modelling has been proposed as a means to control assistive devices by estimating joint torques. Implanted EMG sensors have several advantages over wearable sensors but provide a more localized information on muscle activity, which may impact torque estimates. Here, we tested and compared the use of surface and intramuscular EMG measurements for the estimation of required assistive joint torques using EMG driven modelling. METHODS: Four healthy subjects and three incomplete spinal cord injury (SCI) patients performed walking trials at varying speeds. Motion capture marker trajectories, surface and intramuscular EMG, and ground reaction forces were measured concurrently. Subject-specific musculoskeletal models were developed for all subjects, and inverse dynamics analysis was performed for all individual trials. EMG-driven modelling based joint torque estimates were obtained from surface and intramuscular EMG. RESULTS: The correlation between the experimental and predicted joint torques was similar when using intramuscular or surface EMG as input to the EMG-driven modelling estimator in both healthy individuals and patients. CONCLUSION: We have provided the first comparison of non-invasive and implanted EMG sensors as input signals for torque estimates in healthy individuals and SCI patients. SIGNIFICANCE: Implanted EMG sensors have the potential to be used as a reliable input for assistive exoskeleton joint torque actuation

Voluntary and tremorogenic inputs to motor neuron pools of agonist/antagonist muscles in essential tremor patients

Pathological tremor is an oscillation of body parts at 3–10 Hz, determined by the output of spinal motor neurons (MNs), which receive synaptic inputs from supraspinal centers and muscle afferents. The behavior of spinal MNs during tremor is not well understood, especially in relation to the activation of the multiple muscles involved. Recent studies on patients with essential tremor have shown that antagonist MN pools receive shared input at the tremor frequency. In this study, we investigated the synaptic inputs related to tremor and voluntary movement, and their coordination across antagonist muscles. We analyzed the spike trains of motor units (MUs) identified from high-density surface electromyography from the forearm extensor and flexor muscles in 15 patients with essential tremor during postural tremor. The shared synaptic input was quantified by coherence and phase difference analysis of the spike trains. All pairs of spike trains in each muscle showed coherence peaks at the voluntary drive frequency (1–3 Hz, 0.2 ± 0.2, mean ± SD) and tremor frequency (3–10 Hz, 0.6 ± 0.3) and were synchronized with small phase differences (3.3 ± 25.2° and 3.9 ± 22.0° for the voluntary drive and tremor frequencies, respectively). The coherence between MN spike trains of antagonist muscle groups at the tremor frequency was significantly smaller than intramuscular coherence. We predominantly observed in-phase activation of MUs between agonist/antagonist muscles at the voluntary frequency band (0.6 ± 48.8°) and out-of-phase activation at the tremor frequency band (126.9 ± 75.6°). Thus MNs innervating agonist/antagonist muscles concurrently receive synaptic inputs with different phase shifts in the voluntary and tremor frequency bands

Intramuscular EMG-Driven Musculoskeletal Modelling: Towards Implanted Muscle Interfacing in Spinal Cord Injury Patients

Objective: Surface EMG-driven modelling has been proposed as a means to control assistive devices by estimating joint torques. Implanted EMG sensors have several advantages over wearable sensors but provide a more localized information on muscle activity, which may impact torque estimates. Here, we tested and compared the use of surface and intramuscular EMG measurements for the estimation of required assistive joint torques using EMG driven modelling. Methods: Four healthy subjects and three incomplete spinal cord injury (SCI) patients performed walking trials at varying speeds. Motion capture marker trajectories, surface and intramuscular EMG, and ground reaction forces were measured concurrently. Subject-specific musculoskeletal models were developed for all subjects, and inverse dynamics analysis was performed for all individual trials. EMG-driven modelling based joint torque estimates were obtained from surface and intramuscular EMG. Results: The correlation between the experimental and predicted joint torques was similar when using intramuscular or surface EMG as input to the EMG-driven modelling estimator in both healthy individuals and patients. Conclusion: We have provided the first comparison of non-invasive and implanted EMG sensors as input signals for torque estimates in healthy individuals and SCI patients. Significance: Implanted EMG sensors have the potential to be used as a reliable input for assistive exoskeleton joint torque actuation.The authors would like to thank Enrique Pérez Rizo, Natalia Comino Suárez and María Isabel Sinovas Alonso for their assistance on the experimental and data acquisition procedure

Deep learning for robust decomposition of high-density surface EMG signals

Blind source separation (BSS) algorithms, such as gradient convolution kernel compensation (gCKC), can efficiently and accurately decompose high-density surface electromyography (HD-sEMG) signals into constituent motor unit (MU) action potential trains. Once the separation matrix is blindly estimated on a signal interval, it is also possible to apply the same matrix to subsequent signal segments. Nonetheless, the trained separation matrices are sub-optimal in noisy conditions and require that incoming data undergo computationally expensive whitening. One unexplored alternative is to instead use the paired HD-sEMG signal and BSS output to train a model to predict MU activations within a supervised learning framework. A gated recurrent unit (GRU) network was trained to decompose both simulated and experimental unwhitened HD-sEMG signal using the output of the gCKC algorithm. The results on the experimental data were validated by comparison with the decomposition of concurrently recorded intramuscular EMG signals. The GRU network outperformed gCKC at low signal-to-noise ratios, proving superior performance in generalising to new data. Using 12 seconds of experimental data per recording, the GRU performed similarly to gCKC, at rates of agreement of 92.5% (84.5%-97.5%) and 94.9% (88.8%-100.0%) respectively for GRU and gCKC against matched intramuscular sources

The reliability of methods to estimate the number and size of human motor units and their use with large limb muscles

Purpose: Current methods for estimating muscle motor unit (MU) number provide values which are remarkably similar for muscles of widely differing size, probably because surface electrodes sample from similar and relatively small volumes in each muscle. We have evaluated an alternative means of estimating MU number that takes into account differences in muscle size. Methods: Intramuscular motor unit potentials (MUPs) were recorded and muscle cross-sectional area (CSA) was measured using MRI to provide a motor unit number estimate (iMUNE). This was compared to the traditional MUNE method, using compound muscle action potentials (CMAP) and surface motor unit potentials (sMUPs) recorded using surface electrodes. Data were collected from proximal and distal regions of the vastus lateralis (VL) in young and old men while test–retest reliability was evaluated with VL, tibialis anterior and biceps brachii. Results: MUPs, sMUPs and CMAPs were highly reliable (r = 0.84–0.91). The traditional MUNE, based on surface recordings, did not differ between proximal and distal sites of the VL despite the proximal CSA being twice the distal CSA. iMUNE, however, gave values that differed between young and old and were proportional to the muscle size. Conclusion: When evaluating the contribution that MU loss makes to muscle atrophy, such as in disease or ageing, it is important to have a method such as iMUNE, which takes into account any differences in total muscle size